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ABSTRACT 


This thesis describes an evaluation of the active and 
passive acoustic detection modules of the Naval War College 
Gaming System (NWCGS) that has been installed at the Naval 
War College in Newport, Rhode Island. The specific intent 
of the evaluation is to verify that the model is theoreticaily 
sound. This evaluation compares the NWCGS model to other 
existing acoustic detection models.  Recommendations for 


improvement to the model are also presented. 
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TREE NA NEVER COLLEGE GAMING SYSTEM MODEL 


The detection modules used by the Naval War College 
Gaming System (NWCGS) simulate the acoustic detection proc- 
ess using the generalized expected value sonar equations 
developed during the Second World War. Signal Excess, the 
difference between the existing signal-to-noise ratio and 
that required to achieve a specified level of performance, 
is the fundamental parameter that determines whether or not 
a detection has occurred. In the NWCGS, the value of signal 


excess is determined by the following two equations: 


Ee 


EE SLE MILL DI ND — DT 


when modeling mpassive sonar detections, and 


SE = SL - 2TL + TS - (NL + DT) (2) 


when modeling active sonar detections. 
The parameters referenced by equations (1) anc (2) are 
defined as follows: 
SE - Signal Excess 
SL - The source level of the target as measured one 


yard from the acoustic center of the effective 
sound source 


DI - The directivity index of the hydrophone 
"ED = the one-way transmission loss from the target to 
detector 





NL - The non-target noise in the immediate area of 
the receiver 


F'O be Y RF. strength as measured one yara trom 
the acoustic center of the target 


DT - The detection threshold. This is the signal- 
to-noise ratio (in decibeis) required tu 
achieve a specified probability of detection 
at a specified probability of false alarm 
A more complete description of these parameters is presented 
in Appendix A. The discussion which follows describes the 
methods used in the Naval War College Gaming System to 


determine appropriate values for the parameters of equations 


EV and (2). 


EN METHOD 

The primary method of data retrieval used by tne NWCG 
system is the multi-parameter "look-up" table. A look-up 
table can be described as an array having indices that 
karence its entries. The entries in the look-up table 
describing acoustic transmission loss are transmission loss 
values for different ranges and different enviromrental 
conditions. It is an example of a two-dimensicnal look-up 
table. The first index of the table identifies tle environ- 
mental conditions assigned to the area of the game. The 
second index identifies the distance traveled by sound waves 
from the target to the detector. The advantage of a look-up 
table is its ability to make use of calculations done in 
advance. By using a look-up table, the amount of computations 


done during the course of a game is reduced, and the maximum 


possible game speed is increased. The benefits of a faster 
game are not without some cost. The expense cf using a look- 
up table is the reguirement for additional memory and 
associated hardware used to store the tables. Because a 
look-up table is a discrete approxlriation of a continuous 
process, the accuracy of any signal excess mocel that uti- 
lizes a look-up table is limited by the resolution of the 
tables being used. In the transmission loss look-up table 
example, the look-up table index describing terget-to-sensor 
range has a resolution of 1000 yarcs. This means that 
ranges are rounded to the nearest 1060 yards before entering 
the transmission loss look-up table. Increas:ng the table 
range resolution from 1000 to 500 yards requires twice the 
memory currently required. The presentation that follows 
examines each parameter of equations (1) and (2), and the 


tables used by the NWCGS to determine values Yor each. 


B. TRANSMISSION LOSS 

Transmission loss (described in tne NWCGS documentation 
as propagation loss PL) is determined through a two-parameter 
look-up table. The first parameter determines the environ- 
mental conditions in the area of the game. All players are 
assumed to be contained in a single environmental area, and 
the sound velocity profile is assumed to be the same through- 
out all parts of that area. The second parameter of the table 


is range from target to source (rounded to the nearest 1000 


yards). For a given environmental area, transmission loss 
becomes a function cnly of range. Such a table can adeguately 
account for acoustic spreading losses, however as described 
in Appendix A, absorption and other forms of transmission 
loss depend upon other factors in addition to range. 

l. Direct Path Mode 

When both the target and the searcher are in the 

mixed surface layer, transmission loss is primarily a func- 
tion of range r from target to searcher, target depth DŁ; 


and layer depth D Sound spreads spherically until it fills 


1° 


the layer (at transition range r,), and then spreads cylin- 


E 
drically. Acoustic waves can be scati:ered at the surface, 

or may leak out of the bottom of «she mixed layer. The total 
loss due to surface reflections is equal to the loss per 
surface reflection multiplied by the number of surface 
reflections that occur between tanget and searcher. The 
“ber of reflections occurring is simply the distance between 
target and detector divided by the distance between reflec- 
tions. The complete expression for direct path transmission 


loss (when both target and searcher are within the mixed 


surface layer) is given by the exvression [Ref. 1]: 
TL = 10 logr, t l0 logr + ar + br/r_ (3) 


where: 


r z horizontal range from target to searcher 


= 2 — 
IMSS ED / (D, Di) 


FE 

a. = absorpticn coefficient 

b = loss (dB)/reflection = 0.63*freq* (1.4) sea state) 
= S JS sgrt(D,) 


The transmission loss iook-up table used by the NWCGS model 
describes only losses due to cylindrical spreading. Values 
obtained from such a table can only be used to determine values 
for the first two terms of equation (3). The table does not 
depend upon the environmental area of the ocean (since cylin- 
drical spreading is tne same in all environments); and it 
can be reduced to a one-dimensional table dependent only 
upon the target-to-searcner range. The curve of Figure 1 shows 
how cylindrical spreadiny loss varies as a function of range. 
This plot shows thet a look-up table having a 1000 yard reso- 
lution can be responsipie for table round-off errors that 
exceed 4 dB. In some cases, an error of this magnitude may 
not be acceptable. 

The two transmission loss terms of equation (3) that 
cannot be described by the two-dimensional look-up table 
used in the NWCGS are the losses due to absorption and the 
losses which occur as a result of scattering of acoustic 
energy at the sea surface or "leakage" of energy out of the 
bottom of the layer. Scattering losses are functions of 
range, sea state and signal frequency. A three-dimensional 


look-up table that would be necessary to describe scattering 


10 
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iosses is not an efficient use of available storage, so 
execution time calculation of this loss factor is recommended. 
Some execution time speed can be realized by pre-computing 
the parameter 


E ea) "Sea state (4) 


Tne sea state of a given area remains constant in a given 
area, so b can be precomputed for a number of frequencies 
and stored in a one-dimensional look-up table that is asso- 
ciated with an environmental area. Table l is an example 


cf such a table using a sea state of 4. 


7^ 


2. Absorption Loss 


To account for the additional transmission loss due 
to absorption of sound by the medium, an absorption loss 
factor should be included. The absorption loss corrections 
currently used in the NWCGS model is applied only to the 3 


active model. The correction used is of the form: 
DIM — TE + 0.033 * Free? * range (Kyd) (5) 


This ecuation, cited in Reference 2 is derived from research 
done by M.J. Sheehy and R. Halley [Ref. 3] in 1956. Their 
research attempted to determine the attenuation coefficient 
for low frequency sounds. The experiments involved detonation 


of a nuclear device several hundred miles off the coast of 
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Southern California, and recording the acoustic intensity of 
frequencies between 20 and 200 hertz. The data, collected 
from several locations in the eastern Pacific Ocean, was 
fitted using a least-sguares regression technigue to an 
approximating function described by eguation (5).  Extrapo- 
Lating equation (5) to frequencies above 200 hertz and com- 
paring the predicted results to measured data, their report 
concluded that eguation (5) was a reasonable approximation 
for freguencies up to 60 Kilohertz. 

Subseguent studies, cited by Urick [Ref. 4] have 
developed somewhat different approximating functions for 
the absorption coefficient. The most current of these is 


given by: 


ep ri 2 


2 _—, ee 
SS 


tege 107 dB/Kyd (6) 


A 
H 


Figure 2 is a plot of absorption loss as a function of range. 
The two curves depicted are drawn using two approximating 
functions (eguations 5 and 6). Figures 3 and 4 show how the 
choice of approximating function can affect the total absorp- 
tion loss correction factor for a given range and freguency. 
It is evident from these plots that tne differences between 
the two approximating functions are not large, however as 

the range or freguency increases, the difference between the 
two increases. Assuming that the approximating function 


cited by Urick is based on a larger set of experimental data, 
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it is recommended that eguation (6) replace eguation (5) 
in the NWCGS detection modules. For the sake of consistency, 
it is further recommended that the absorpticn loss correc- 
tion be applied to acoustic signals in the passive acoustic 
detection modules as well as the active acoustic detection 
modules. In the case of active detections, the frequency 
"£f" of equation (6) is the ping frequency of the sonar trans- 
ducer. In the case of passive detections, "£" should be the 
center frequency of the bandwidth in which tne searcher is 
looking. Since equation (6) is a function of frequency 
alone, a one-dimensional absorption coefficient iook-up table 
can be constructed. Table 2, derived from equation (6), 
provides data that may be useful in the construccion of such 
a look-up table. 
See convergence Zone Mode 

If convergence zone transmission paths are possible, 
the NWCG system is able to adjust the transmission loss to 
reflect this type of contact. Whenever the searcher is in 
the target's convergence zone area, the NWCG system multi- 
plies the unadjusted transmission loss by a convergence zone 
correction factor. In order to reduce the unadjusted trans- 
mission loss to reflect a convergence zone gain, the multi- 
plicative factor must be a real number between zero and one. 
This method is significantly different than the method used 
by Coppens [Ref. 5] to describe this gain (see Appendix A). 


Since this method of adding rather than multiplying the 
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correction factor has become an accepted standard, it is 
recommended that tne NWCG system be modified to take advan- 
tage of the availability of existiny data bases. 
4. Bottom Bounce Mode 

Transmission loss is adjusted in the NWCG system 
whenever the environment supports bottom bounce. A bottom 
bounce flag is set if the range from target to searcher is 
within a bottom bounce region, and all transmission loss 
calculations are multiplied bv a bottom bounce correction 
factor. The areas where bottom bounce detections are con- 
sidered likely are those between the maximum direct path 
range and the minimum convergence zone range. While bottom 
bounce transmission paths can extend into both direct path 
and convergence zone areas, the trensmission loss for bottom 
bounce is generally larger thar that of the other two paths. 
Consequently, any detection thet occurs in one of the other 
transmission areas will most likely not be due to bottom 
bounce. The current literature discussed in Appendix A 
describes the following procedure to adjust the transmission 
loss: 

First, increase the effective range from target to 


searcher by the eguation: 


R = R/(cos 9) (7) 


where 0 is the angle made between the incident sound wave 


and the bottom. The increase in range is due to the indirect 
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transmission path. Second, increase the total transmission 
loss by adding a constant factor BL due to absorption and 
scattering of the signal by; the ocean bottom. 

The spreading ioss associated with bottom bounce 
transmission paths is spherical. Since the NWCGS spreading 
loss table describes cylindrical spreading, spreading loss 
values obtained from this table must be doubled to reflect 
. that difference. The loss associated with a bottom reflec- 
tion is an environmental parameter that should be included 
in the table that describes the environment. 

An important aspect of bottom bounce often overlooked 
is its value in cases of layer depth interdiction. When the 
searcher is above the mixed layer, and the target is below 
that layer (or vice versa) direct path ranges are severely 
limited due to the reiracting effects of the layer. Bottom 
bounce signals have much Figher angles of incidence, and are 
therefore less affected by refraction. If layer depth inter- 
diction occurs in bottom bounce areas, detections should be 
based strictly on bottom bounce transmission paths. 

Sen conclusion 

Figure 5 shows typical transmission loss profile for 
the three transmission modes (direct path, bottom bounce, 
and convergence zone). Using the equations of the NWCG sys- 
tem, a Similar transmission loss profile has been computed 
and shown in Figure 6. The two figures are significantly 
different, especially in terms of bottom bounce transmission 


loss. A more extensive transmission loss look-up table is 
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reguired to adeguately describe spreading, absorption, trans- 
mission path, etc. This table could eliminate the need for 
bottom bounce and convergence zone correction factors, and 
could be used to describe environmental phenomena that are 
difficult to describe analytically. The NWCGS transmission 
loss look-up table could be modified to include frequency, 

in addition to environmental area and range, as an entering 
argument to the table.  Additionally, the resolution of the 
table (currently 1000 yards) could be increased. 

If the current structure of the transmission loss 
look-up table cannot be altered, the following recommenda- 
tlons are made: 

l. Use the propagation loss table currently in the NWCG 
system to describe transmission loss due to cylindri- 
cal sreading only. Spreading losses are a function 
of range only, and the difference between cylindrical 
and spherical spreading losses is a multiplicative 
factor of two. All transmission paths are subject 
to spreading losses, and therefore all transmission 
paths will reference this table. 

2. For direct path transmissions, transmission loss should 


be described by the equation: 


TL = l0logr, + 10 logr + ar + br/r, (3) 


where: 


24 


s 2 
te = 115 sgrt (D,/(D, -D ) 


E tl 
Eege e 918 sgrt(D,) 
a = absorption loss/yard 
b = loss (dB) /reflection 


pee SCE), freg) 


The value r. is fixed by the environment selected. 
Consequently, it is possible to compute the value of 

10 log (r,) beforehand, and include it in the parameter 
list of the environmental area. The value of 10 log (r) 
is obtained directly from the propagation loss look-up 
table. Absorption loss per yard "a" can be determined 
by a call to another look-up table (see Appendix A). 
The loss per bounce parameter "b" is a function of both 
sea state and frequency; and a two-dimensional look-up 
table similar to Table 2 can be constructed and used 
to determine "b". 


For bottom bounce transmission paths, the transmission 


loss should be reflected by the equation: 
TL = 20 log (r/cos 8) + ar/cos 8 + BB (99 


The range from target to searcher is divided by cos 9 
to account for the increased transmission path. The 


adjusted range is an entering argument for the propagation 
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loss look-up table. Since bottom bounce spreading 

is spherical rather than cylindrical, the value 
obtained from the table should be doubled. Absorption 
coefficient "a" is the same as previously described, 
and the bottom bounce loss BB is a parameter determined 
by the environmental area. 

4. In areas wihtout bottom bounce, if the detector is 
above the mixed layer and the target is below (or vice 
versa), it is reasonable to assume no detections are 
possible. Under the same circumstances in a bottom 
bounce area, detections via bottom bounce transmission 
paths are possible, and should be investigated. This 
would require the modification of the layer depth 
interdiction test. If the detector is above the laver 
and the target is below (or vice versa) and the bottom 
bounce available flag is not set, then no detection 
Can occur. If the bottom bounce available flag is 
set, transmission loss is calculated for bottom bounce 


mode. 


o TRRESHOLD CROSSING MODELS 

A major use of the passive or active sonar equation is 
to determine mean signal excess. Mean signal excess is the 
average difference (in decibels) between the measured signal- 
to-noise ratio and the signal-to-noise ratio which will yield 
a 50 percent probability of detection for a specified proba- 


bility of false alarm. The NWCGS utilizes a signal excess 
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threshold crossing model to account for the randomness in 
the components of the sonar eguation, and to determine the 
probability of detection. This model defines tne signal 
excess at a time t to be a random variable Kae LL and a 
detection to be the event (X..(t) > 0}. The signal excess 


SE 


at time t can be expressed as: 


Xop (€) = SE(t) + X(t) (11) 


where SE(t) is the mean signal excess, and X(t) is a random 
variable that determines the variability of the signal 
excess at time t. Values for X(t) are determined by drawing 
from a normal distribution having a mean of zero, and a 
varlance that is determined by the variability of the com- 
ponents of the sonar eguation. The probability that a 


detection will occur can be written as follows: 


P(X--(t) > 0} = PiX(t) > -SE(t)) (12) 


Since X(t) is normally distributed, this implies that 


PiK.p(t) > 0) = {SE(t)/o} (150 
The mean signal excess is used as an entering argument for 
an inverse normal look-up table that determines the proba- 


Pretty of detection. 





Two modes are available for converting the detection 
probability into a detection decision, the deterministic 
mode and the probabilistic mode. In the deterministic mode, 
the probability of detection is compared to a probability 
of detection threshold that has been set by the umpire during 
pne initialization of the game. If the computed probability 
of detection exceeds this predetermined threshold, a detection 
is recorded. The alternative to the deterministic mode is 
the probabilistic mode. In this case, the probability of 
detection is computed, and the detection process is simu- 
lated by a draw from a Uniform (0,1) distribution. a detec- 
tion occurs whenever the random variable drawn is less than 
the computed probability of detection. 

When using the deterministic mode with a giver standari 
deviation o, the probability of detection threshold ir.itlalized 
by the game director defines a signal excess to signa ratio 


threshold (SE/c) by: 


TA 


$(SE/0) (PA) my TE) 


TH 
If the probability of detection threshold is initialized to 


0.5, then (SE/o) = 0 and a detection occurs anytime the 


TH 
mean signal excess is equal to or greater than zero. In the 
probabilistic mode, the value ${SE/o} is compared to a random 


draw. If o > 0, then for any non-negative value of SE, 


GEE) > I, and detections occur with probability 1. For 
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any negative value of SE, ${SE/o} > 0, and detections occur 
with probability 0. Consequently, in the limiting case where 
o = 0, detections occur whenever the mean signal excess is 
equal to or greater than zero, and in this limit, the 
probabilistic mode is equivalent to the deterministic mode 
with the probability of detection threshold set to 0.5. This 
suggests that the deterministic mode can be replaced by this 
limiting form of the probabilistic mode unless there is some 
MeLLity in having (Pd) my set to a value diiferent from 0.5. 
The probabilistic mode corresponds to a complete inde- 

pendence model. In this mode, detection on a looh is inde- 
pendent of detection on any other look. The cumu..ative 
probability of detecting a target in N looks is given by: 

N 


P(detecting in N looks} = 1- mI (l-Pd(i); 
isl 


1 (15) 


where PD(i) is the detection probability on the ed look. 
For any non-zero probability of detection Pd, the cumulative 
probability of detection approaches one as tne number of 
looks increase. Figure 7 shows the cumulative probability 
of detection for a complete independence model as a function 
of the number of looks, where the mean signai excess is 
constant and the single look probability of detection is 
OI, 


This figure also indicates the sensitivity of this model 


to the time between looks. It shows that after 22 looks, 
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the cumulative probability of detection exceeds 90 percent. 
It can also be shown that the average number of looks to a 
detection is 10. For this mcdel to be valid, the time be- 
tween looks must be long enough to have the fee pendence 
assumption hold. In the case of passive acoustic detections, 
the time between looks must also be long enough to satisfy 
the integration time requirements for detection threshold. 
In the NWCGS, the time between looks is a director adjusta- 
ble parameter initialized to 20 minutes. It is important 
that the game director understand how changing this value 
affects the cumulative pbrobabili-y of detection. The dis- 
cussion which follows descripes some alternative models. 
ie (AsO) Jump Model 

A more complicated model for the random component of 
signal excess than the one used in the NWCGS was introduced 
by J.D. Kettelle in 1960. it is referred to as the (1-0) 
jump model. In this model, the initial value of X(t) is 
determined in the same manner as that used in the NWCGS 
model; that is by a draw from a normal distribution having a 
mean of zero and a standard deviation of o. On subsequent 
looks, the value of X(t) remains the same until a jump occurs. 
The time of the jump is determined by a draw from an exponen- 
tial distribution having a mean of A. Each time a jump occurs, 
a new value for X(t) and a time to the next jump are deter- 
mined. The choice of lamda (known as the inverse relaxation 


constant) determines the model's level of dependence. As 
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A > ©, the mean time between jumps approaches zero, and the 
model apprcaches the independence model. As A > 0, the time 
between jumps approaches o, and this limiting case corresponds 
to a complete dependence model. Complete dependence refers 
to the fact that if a detection does not occur on the detection 
opportunity, or "lock" with the largest mean signal excess, 
then detection will not occur on any other look. Therefore, 
the cumulative probability of detection for a series of looks 
is equal to the probability of detection for the look with 
the largest mean signal excess. Figure 8 shows a number of 
possible cumulative probability of detection curves available 
for a given mean nrcbab:.lity of detection (0.1) and different 
values for A. 

2. Gauss-Markov Models 


An alternative to the (A-c) jump model is the Gauss- 


Markov model. This model is Gaussian in nature because the 
E distribution of a set of looks [X(1),X(2),...,X(n)] is 
we variate normal distribution. For any X(i), its 


marginal distribution is normal with a mean of zero, and a 


variance of EE It is Markovian, implying that given values 
mere K(0O), X{(l;, ..., X(t-l), the value of X(t) depends only 
upon X(t-1). This conditional dependence of X(t) on X(t-l) 


is normally distributed with a mean of o*X(t-1), and a variance 


2 


of See ). The parameter p, known as the autocorrelation 


coefficient, is defined by the equation: 


n= e CES) 
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MA NPABENCEST ACT Equation (15) corresponds to the inverse 
relaxation coefficient of the (A-c) model. As A > 0, the 
parue cf o > l and, as with the (A-o) model, the Gauss-Markov 
model approaches a complete dependence model. Similarly, 

as A > o, the value of p > 0 and the Gauss-Markov model 
approaches the independence model. The degree of dependence 
ie determined by choice of A. 

Although the Gauss-Markov model is difficult to deal 
with aralytically, a computer model that uses the conditional 
EE LrFYIEbutron of X(t) given X(t-l) is not difficult to develop. 
Cumulative probability of detection curves derived through 
simulation using various values of lamda are presented in 
Figure 9. 

A Out of N Models 

The k out of N detection model that is described here 
is based on independence model, however, each look results 
in either a success or failure rather than a detection or a 
ncn-dezection. For a detection to occur, k successes must 
have occurred within the last N looks. Such models provide 
a means of describing operator characteristics. Figure 10 
shows how the k out of N model compares with the dependence 
and independence models of the NWCGS. The k out of N model 
offers the advantage of allowing one to set the time between 
looks to any specified time interval and then, by choice of 
k and N, of fitting the model to an appropriate cumulative 


probability of detection curve. Most detection processes 
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reguire an initial integration time before a detection can 
be called, and this model can easily accommodate such a 
phenomenon. A k out of N model could also have been based 
on either the (A-o) jump model or the Gauss-Markov model. 
An example of a more general k out of N model can be 
found in the APSURF Surface Ship Engagement Model [Ref. 6]. 
In that model, active sonars were modelled as follows: 
l. A 6 out of 28 model determined detections if the 
operator was unalerted. 
2. A 5 out of 20 model determined detections if the 
operator was alerted. 
3. A 2 out of l0 model was used after a detection to 
determine if contact was still being maintained. 
This model is easily implemented and would require very 
little in the way of software modification to the existing 
NWCGS model. 
4. Summary 
The alternative models described above are a sample 
of the models that have been used to describe the detection 
process. Although none of the models are without shortcomingz, 
each of them has desirable properties not found in the current 
NWCGS models. In particular, they provide a range of serial 
correlation properties. This is appealing in modeling a 
detection process where the probability of detection time t 


may be influenced by information obtained at earlier times. 
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APPENDIX A 


The most effective mechanism for locating targets under 
the ocean surface is the detection and evaluation of acous- 
tic information specific to a target of interest. With 
this mechanism, detection may occur as a result of passively 
listening to the sounds generated by the Faro šuje propulsion 
equipment and auxiliaries, or may involve actively radiating 
acoustic energy in the suspected direction of the target anč 
listening for the echo which would reflect off the target 
and return to the searcher. Passive sonar systems otfer the 
advantage of being able to detect noisy targets at extremely 
long ranges (often in excess of one hundred miles) without 
alerting the target to the starch. Although active sonar sys- 


tems, or echo ranging systems, alert the target to what the 


fu 


searcher is about, they can be most effective against 
quiet, slow-moving target close to the searcher. The system 
of choice is determined by the tactical situation and the 
environmental conditions under consideration. 

The passive acoustic modules used in the NWCGS are based 
on the generalized expected value sonar equations first for- 
mulated during the Second World War. Signal Excess, or the 
difference between the existing signal-to-noise ratio and 
that required to achieve a specified level of performance 


[Ref. 7] is the fundamental parameter that determines whether 


or not a detection has occurred in these modules. The required 
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signal-to-noise ratio (in decibels) is cailed the detection 
threshold. It is defined to be the signal-to-noise ratio 
that will give a specified probability of detection for a 
specified false alarm probability in a forced choice decision 
Situation. If the specified probability of detection is equal 
to 0.5, detection threshold is usually referred to as recog- 
m@eeron differential. 

After determining the appropriate detection threshold, 
the detection model determines the current signal excess 


level using the formula 

SE = SL - IL = NE + DI = DT (21) 
for modeling passive detections, and 

SE = SL - 2TL + TS - NL + DI - DT (2,2) 


for modeling active detections, where the parameters of the 


equations are defined as: 


SE - Signal Excess 

SL - the source level of the target as measured 
one yard from the acoustic center of the 
source 

DI — The directivity index of the receiver 

TL - The one-way transmission loss from the target 


to the source 


NL - The non-target noise existing in the area of 
the receiver 
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IS - The target strength as measurec one yard from 
the acoustic center of the target 


DT - The detection threshold, or the minimum signal- 


to-noise ratio (in decibels} that is required 
to call a detection 


A more detailed discussion of these parameters follows. 


me SOURCE LEVEL 

The value of the acoustic energy radiated iato the medium 
by the searcher is termed the projector sonrce level. For 
the sake of standardization, the source ievel is taken to be 
the acoustic power as extrapolated bacx to vne yard from 
the acoustic center of the sound source. The interference 
patterns and other physical characteristics of acoustic 
energy close to its source produce measurements which are 
OE little ine, SO in practice the acoustic measurements 
are made from large distances, and the results can be 
reliably extrapolated to give reasonable ssurce levels for 
all meaningful ranges. In passive acoustic listening, the 
Signals of interest are those which are generated by the 


target alone. 


5, DIRECTIVITY INDEX 

Total acoustic power arriving at the receiving unit (be 
it a single hydrophone or array of hydrophones) consists of 
both signal power and noise power. Signal power is generally 
assumed to arrive at the receiver from a signal direction. 


Conversely, noise power is generally assumed to arrive at 
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the receiver uniformly from ail directions. A receiving unit 
that only accepts acoustic power from specific directions 
usually has a higher signal-to-noise ratio than a recelving 
unit that accepts acoustic power omni-directionally. The 
signal-to-noise ratio gain due to the directivity of the re- 
ceiving unit is termed the directivity index DI of the 


recelver. 


C. TRANSMISSION LOSS 

As with any other form of radiated energy, acoustic 
Signals are subjected to a number of paysical effects which 
tend to distort and dissipate the signal. In the acoustic 
sonar models, these effects are collectively described by 
the parameter transmission loss (TL). The primary constitu- 

e 
ents of transmission loss are spreading, and absorption. 

An omni-directional sound sou:zce in a lossless infinite 
medium will radiate the same acoustic power over an ever- 
increasing area. Acoustic intensity, the quantity being 
measured by a detection system, is equal to power "P" divided 


by the total area "A" over whicn the intensity is distributed: 
P = I*A (A5) 


If acoustic energy is subject only to spherical spreading in 


a lossless medium, the transmission loss can be described by: 


TL. OO EE) (A6) 
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where r is the radius from the center of the sound source. 
In this case, intensity is inversely proportional to the 
sguare of the distance between detector and the center of 
the sound source. 

When physical or environmental boundaries serve to inhibit 
spherical spreading by restricting sound between two parallel 
planes (the ocean floor and surface as an example), the power 
is radiated more like a cylinder of constant height "h" and 
increasing radius "r". The area of this cylindrical acoustic 


wave would be: 


"X em 2*n*p^*h (A7) 


Intensity is inversely >rcportional to the first power of 
range making the spreading loss component of transmission 
loss equal to 


TL. = 10 % LOG (5) (A8) 


Initially, the intensity of sound sources spreads spherically. 
At some distance r, away from the sound source, the acoustic 
intensity fills the medium between its bounds and spreading 
losses become more cylindrical than spherical. If the 
searcher is further away than this transition range, the 
transmission loss equation can be described by 


ier. 10 * (LOG LR, 7210.27 LOG (r) (A9) 
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The most common instance of this type phenomena both the 
target and the searcher are in the surface layer where the 
sound velocity profile is such that a lower acoustic "boun- 
dary" is formed. A similar, though less common example is 
the SOFAR, or deep sound channel formed when the target and 
searcher detection system are both located in a channel 
where sound is trapped from above and below by acoustic 
"boundaries" [Ref. 8]. In either event, the threshold range 
R, is related to the height of the channel D 


t 
target below the top of the channel D 


1 and depth of 


t PY 


R. = 1.094 * sqrt (D{/(D,-D,)) (ALO) 


e E 


where all quantities are measured in yards. 
1. Convergence Zone 

The nature of the oceans are such that given suffi- 
cient depth, all acoustic signals will be refracted back 
toward the surface. In this refraction process, the sound 
waves tend to ccnverge, thereby reducing the losses due to 
Spreading. These convergence rings, or zones as they are 
called, have transmission losses at least 10 dB less than 
they would be if the transmission were only spherical spread- 
ing. To account for this, transmission losses are calculated, 
and a convergence zone gain parameter CZ is subtracted from 


the expected transmission loss. 
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2. Bottom Bounce 

The mud-sand composition of much of the ocean floor 
tends to absorb rather than reflect acoustic signals. If 
the bottom consists more of rocks and sand than mud, the 
reflective properties of those materials will allow trans- 
mission of sound at greater than anticipated ranges. To 
account for this effect, a bottom bounce adjustment factor 
BB can ce subtracted from the transmission loss calculation. 

3. Absocptlon Losses 

Afte:ic spreading losses, the second greatest cause of 
intensity dissipation is absorption. This is caused by the 
conversion of acoustic energy into another form of energy. 
This type loss is usually attributed to shear viscosity, a 
friction-like conversion of acoustic energy into heat, or to 
energy involved in the dissociation and reassociation of 
ions in the ocean. The shear viscosity is associated with 
all frecuencies of sound, while the dissociation reassocia- 
tion ohénomena predominate at frequencies below 100 kHz. 
The intensity if absorbed exponentially with distance travelled 
through the absorbing medium and transmission losses asso- 


ciated with absorption can be described by: 
IL. = ar (All) 


where a is the absorption coefficient of the medium in units 
of decibels per kiloyard. A number of studies have been 


done in an effort to define an appropriate value for alpha. 


44 





it has been shown that the [Ref. 9] absorption coefficient 
is primarily a function of frequency and can be approximated 


using the equation: 


2 


a= 1.86 x 10 — + 2.86 x 10 =— dB/Kyd (A12) 
2 2 F 
Ea T f T 


where f is the frequency of the sound, S is the salinity 
(in parts per thousand) of the sea, and Fm is a relaxation 


freguency related to temperature T by the eguation: 


C=|FSZWACT Z 737P1) 


FP, = 21.9 * 10! (A13) 


di 


Additional causes of transmission losses are absorption of 
sourd into the ocean floor, scattering of sound energy at | 
the surface, leakage of sound energy out of a sound channel, 
changes in the absorption coefficient due to depth and tem- 
perature changes, etc. A complete accounting for all the 
elements which comprise transmission losses would be compu- 
tationally impractical and unnecessary for the model being 
discussed. As a minimum however, losses due to spreading 
and absorption should be considered. 

Transmission loss models are primarily a function of 
frequency of the sound source, and distance from the source 
to the target. As such, the one-way transmission loss from 
target to detector will be the same for active and passive 


sonars (assuming both systems are using the same frequency). 
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Active sonar signals must travel from source to target and 
back, so the transmission loss in that circumstance is 


assumed to be exactly twice that of a passive sonar signal. 


D: NOISE LEVEL 
Any sound not originating from the target of interest 
is considered to be noise in the sonar model. Noise levels 
are generally categorized as either ambient noise (AN) or 
searcher self-noise (SN). 
Ambient noise is the inherent noise of the ocean itself. 
It is a collection of acoustic signals emanating from numer- 
ous sound sources randomly arriving at the location of 
the searcher's hydrophone. It is primarily a result of the 
collowing factors: 
l. Environmental effects--tidal activity, rain, wind, 
seismic activity, turbulence, sea state, etc. 
2. Biologics--whales, dolphins, snapping shrimp, 
croakers, etc. 
3. Man-made effects--shipping intensity, petroleum 
exploration, deep water mining, etc. 
The diversity of sound sources involved in the makeup of 
ambient noise means that the entire acoustic spectrum is 
affected. It also implies that ambient noise is likely to 
be one of the most variable factors in the sonar equation. 
Ambient noise is usually modeled as a normal random varia- 
ble whose mean is dependent upon the shipping intensity in 


the area of the searcher and target [Ref. 10]. 
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Another major contributor to the undesirable noise at 
the detector is the self noise of the searcher.  Auxiliary 
eguipment is always being operated, and propulsion eguip- 
ment is nearly always in use. As the speed of the searcher 
increases, so does the self noise problem associated with 
the propulsion equipment. At higher speed, the turbulence 
created by water flowing around the hydrophones becomes an 
additional source of self noise.  Hydrophones are generally 
placed well forward of machinery spaces and have designed 
directivity away from any self noise sources, however total 
elimination of self noise is impossible. Self noise is 
usually modeled as a function of the speed of the searcher. 
ERY" TARGET STRENGTH 

When a target is ensonified by an dd" detection sys- 
tem, a portion of the incoming power is reflected back to 
the searcher. If the target were a perfectly reflecting 
spherical body, the sound energy would be reflected uni- 
formly about the target, spreading spherically. In practice, 
targets are neither spherical, nor perfectly reflective. 
Target strength (TS) is the parameter that accounts for both 
absorption of energy by the target, and its reflectivity in 
the direction of the searcher. Target strength is strongly 
influenced by the relative angle of the target and searcher 
(aspect angle). The variance of target strength with aspect 
differs widely among targets, but as a first order approxi- 


mation, Urick [Ref. ll] proposed an average submarine 


47 





target strength to aspect relationships which can be des- 


cribed by the eguation: 


TS = TS "(16.17-2.98 cos (28)-3.083 cos (68) (A14) 


where TS, is the maximum value of target strength, and B 
is the aspect angle between target and searcher. Figure ll 
shows how this eguation relates target strength to aspect 


angle. 


RR DETECTION THRESHOLD 
Target information available to the searcher is processed 
in some manner to arrive at a detection decision. The signal 


= 


excess models used in the NWCGS assume that detections occur 
whenever the signal-to-noise mo (in decibels) exceeds sore 
predetermined detection threshold DT. The choice of DT 
determines the probability of making a decision error. In 

a situation where a system is forced to make a decision after 
a single look, a detection will be called if the input exceeds 
a threshold. When no target is present, the input is due 

to noise alone. In this event, if the input level exceeds 
the threshold, a detection will be erroneously calied, and a 
decision error known as false alarm occurs. When a target 

is present, the input level is both signal and noise. In 


this event, if the combined signal and noise input level 


exceeds the threshold, the system will correctly decide that 
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a target exists.  Specifying the acceptable probacility of 
false alarm determines a threshold. Having established a 
threshold, a specified probability of detection determines 
a required signal-to-noise ratio. This signal-to-noise 


ratio (in decibels) is the detection threshold DT. 
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